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ABSTRACT 

The internal energy available in vibrationally excited H 2 molecules can be used to overcome or 
diminish the activation barrier of various chemical reactions of interest for molecular astrophysics. 
In this article we investigate in detail the impact on the chemical composition of interstellar clouds 
of the reactions of vibrationally excited H 2 with C + , He + , O, OH, and CN, based on the available 
chemical kinetics data. It is found that the reaction of H 2 (v > 0) and C + has a profound impact on 
the abundances of some molecules, especially CH + , which is a direct product and is readily formed 
in astronomical regions with fractional abundances of vibrationally excited H 2 , relative to ground 
state H 2 , in excess of ~ 10~ 6 , independently of whether the gas is hot or not. The effects of these 
reactions on the chemical composition of the diffuse clouds C Oph and HD 34078, the dense PDR 
Orion Bar, the planetary nebula NGC 7027, and the circumstellar disk around the B9 star HD 176386 
are investigated through PDR models. We find that formation of CH + is especially favored in dense 
and highly FUV illuminated regions such as the Orion Bar and the planetary nebula NGC 7027, where 
column densities in excess of 10 13 cm -2 are predicted. In diffuse clouds, however, this mechanism is 
found to be not efficient enough to form CH + with a column density close to the values derived from 
astronomical observations. 

Subject headings: astrochemistry — molecular processes — ISM: molecules, photon-dominated-region 
(PDR) 



1. INTRODUCTION 

As concerns the chemical modelling of interstellar 
clouds, state-to-state chemistry, i.e. that which con- 
siders the quantum states of reactants and/or products 
in a given chemical reaction, has been traditionally ig- 
nored, mainly because the required chemical kinetics 
data are generally not available and because it is ex- 
pected that it would not affect in a drastic way the 
predicted chemical abundances. Some attention how- 
ever has been given to the different reactivity of or- 
tho and para states of molecular hydrogen and its effect 
on various aspects of interstellar chemistry such as the 
formation of ammonia or the ortho-to- para ratio of c- 
CsHg in cold dense clouds (|Le Bourlotlll99H iPark et all 
2006). The enhancement in the reactivity of molecu- 
lar hydrogen when it is in a vibrationally excited state 
has also been considered in the past, mainly in the con- 
text of trying to explain the long standing problem of 
CH + formation in diffuse cloud s dStecher fc Williams! 
H97l iFreeman fc Williams! [1981: IGarrod et alT T2003L 
but also in the study of the che mistry of dense pho- 
ton dominated regions (PDRs) dTielens fc Hollenbachl 
119851: iSternberg fePalgarnol 11995ft and shocked region s 
(|Wagner fc Graffl 119871 : IHollenbach fc McKeel Il98flh . 
These studies assumed different reactivity enhancements 
for reactions of vibrationally excited H 2 with species such 
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as C + , O, or C, and found that whereas in the case of dif- 
fuse clouds a too large fraction of H 2 v > is required to 
explain the formation of CH + with the observed abun- 
dance, in the case of dense PDRs there are important 
effects for the chemistry. 

In this paper we look again at the chemistry of vi- 
brationally excited H 2 in the light of some experimen- 
tal and theoretical studies carried out in the last years, 
which have been essentially ignored by the astrochemical 
community, and investigate the implications for various 
astronomical regions where a relatively large fraction of 
vibrationally excited H 2 has been observed or is expected 
to be present. 

2. THE CHEMISTRY OF VIBRATIONALLY EXCITED H 2 
2.1. Basics 

The main feature that characterizes the chemistry of 
H 2 in excited vibrational states -as compared to that 
of H 2 in the ground vibrational state- is that the in- 
ternal energy of the excited hydrogen molecule can be 
used to overcome or diminish activation barriers which 
are present when H 2 is in its ground vibrational state 4 . 
Several experimental and theoretical studies have demon- 
strated that such behavior occurs for some reactions of 

4 In regard to the chemistry of vibrationally excited H2 it is 
worth to comment on the difference between a state-specific and a 
thermal rate constant. The thermal rate constant, usually provided 
by experimental and theoretical chemists and used in interstellar 
chemical models, is the thermal average of the individual state- 
specific values, i.e. the sum over all H2 states of the state-specific 
rate constant multiplied by the thermal fractional population of 
that H2 state. It turns out that since the vibrationally excited 
levels of H2 are very high in energy and thus poorly populated 
in thermal conditions, except at high temperatures (> 1000 K), 
we expect the thermal rate constant to be essentially equal to the 
state-specific H2 v = rate constant, so that hereafter we use both 
terms interchangeably. 
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TABLE 1 

Thermal and state-specific rate constants for chemical reactions of H 2 relevant for astrophysics 
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H 2 + H 
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1.52X10" 11 (T/300) 1 - 33 exp(-902/T) 
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-13360 
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HON + H 


1.17X10" 12 (T/300)^- 31 exp(-1188/T) 
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H 2 (u 
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CN ->• HCN + H 


9.65X10- 12 (T/300) 1 - 04 exp(-1397/T) 


200-1000 
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Re fere nces. — (l) I Hierlet alj 1 119971 ); (2) based on IGerlich et all J1987D ; (3) UMIST Database for Astro chemistry <|Woodall et all 
[2007D; (4) [Jonesjet_ajj y.986D ; (5) NI S T Chem i cal K ineti cs Database (http ://kinetics.nist.gov/kinetics/); (6) Suit an ov &c Balakrishnan 
1120051) ; (7> IZellner & Steinerfl h9811hlTruon3 J 19951) ; (8) |Zhu et all Jl998T) . 

Note. — When no specific state of H2 is indicated the expression for k corresponds to the thermal rate constant. The reaction enthalpies 
AH® (OK) ha ve been computed from the formation enthalpies at K of each species, taken from N IST- J AN AF Thermo chemical Tables 
flChase||199"g|) , while for the H 2 state-s pecific reactions the level energies of H2 have been taken from [Dabrows 
a AT is the temperature range over which the rate constant has been studied. 

b Ej is the energy of each rotational level of H2 (v — 0, j) relative to the ground state V — 0, j — {E\ — 170.5 K, E2 — 509.9 K, E3 — 1015.2 
K, £ 4 = 1681.7 K, _E 5 =2503.9 K, _E 6 =3474.4 K, and _E 7 =4586.4 K : I D abrowskill 1 9841) . 

c In the temperature range 500-4000 K the rate constants for the reactions between O and H2 in the v — 1 and v — 2 states can be better 
fit by the expressions 7.79X10" 12 (T/300) 1 20 exp(-2444/T) and 2.39X10 -11 (T/300) 87 exp(-1325/T) respectively. 



H2 which are of relevance for astrophysics, the kinetic 
data of which are detailed in Tableland Fig. [T] 
The hrst reaction that draw our attention is 



The second reaction in which we focus is 
H 2 + He+ -> He + H + H+ 



Ho 



CH 



H 



(1) 



which is well known to possess an endothermicity of 0.37 
eV (~ 4300 K). It has long been thought that such en- 
dothermicity could be overcome by vibrationally excited 
H2 (the v = 1 levels have energies > 0.5 eV above the 
ground state), hypothe sis that gained exp erimental sup- 
port with the study of Uones et al.l (1 19861) and the more 
recent measurements of iHierl et al.l (|1997l ). These lat- 
ter authors were able to measure the rate constant of 
H2(w = 1) + C + over the 800-1300 K temperature range 
and found it to be in the range 1-2 xl0~ 9 cm 3 s _1 , i.e. 
essentially equal to the Langevin collision rate (1.6 xlO -9 
cm 3 s _1 ), and to show a slight negative temperature de- 
pendence. This is the typical behavior of exothermic 
ion-non polar neutral reactions which have a temper- 
ature independent rate constant close to the collision 
limit, and indicates that it is quite safe to adopt also 
the collision rate at low temperatures for the reaction 
of H2(w = 1) and C + . In connection with this we note 
th at rate coefficients fo r reaction ([I]) have been computed 
by IGerlich et al.l (|1987D for selected initial rotational lev- 
els (.7 = 0, 7) of H2 in the ground vibrational state. This 
work was combined with guided-beam measurements and 
a good agreement between the experiment al results and 
the st atistical calculations was observed. IGerlich et al.l 
(1987) found that the state-specific rate constants may 
be well approximated by Arrhenius expressions (see Ta- 
ble [T] and Fig. QJ in which as j increases the exponential 
factor decreases by an amount equal to the energy of the 
j level. All this indicates that the energy of excited H2 
in reaction (JIJ is effectively used to diminish or overcome 
the endothermicity of the reaction, and that for H2 states 
with v > the reaction proceeds at the collision rate. 



(2) 

which is strongly exothermic but shows an extremely 
small rate constant, some 10~ 14 cm 3 s^ 1 below 300 K 
l|Schauer et al.|[i"989l) , which in dicates the p r esence of an 
important activation barrier. Uones et"aLl (|1986f) have 
studied the kinetics of this reaction experimentally at 
room temperature and found that when H2 is in a vi- 
brational state v > 1 the reaction proceeds at about the 
collision limit, 1.8 xl0~ 9 cm 3 s _1 , value that we adopt 
hereafter. 
Another interesting reaction is 



H 2 + O -+ OH + H 



(3) 



which is known to be endothermic by 0.08 eV (~ 900 K) 
and to possess an even higher activation barrier of about 
0.4 eV (~ 4800 K). Recently ISultanov fc Balakrishnan! 
(2005) have calculated the rate constants for reaction 
(3]) with H2 in various vibrational states (v = — 3) and 
have found that although the rate constant is noticeably 
enhanced as the vibrational state of H2 is increased, the 
reaction does not proceed at the collision rate limit for 
H2 (v > 0) but an activation barrier persists even when 
H 2 is in the v = 3 state (see Fig. Q}. 
The exothermic reaction 



H 2 + OH^H 2 + H 



(4) 



possesses an activation barrier which can be partially 
overcome, when H2 is in first excited vi brational state 
(|Zellner fc Steinertl Il98l ITruonel fl995L see Fig. Q). 
Something similar occurs for the exothermic reaction 



H 2 + CN -+ HCN + H 



(-5) 



the rate constant of which experiences a moderate 
enhancement w hen H2 is in the vibrational state v = 1 
(|Zhu et al.lll998l see Fig. [J). 
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Fig. 1. — Rate constants as a function of the reciprocal of tem- 
perature for thermal and state-specific reactions of H2 with C~*~, 
He+, O, OH, and CN. 

The studies carried out on the reactions commented 
above leave us the following general picture. The reac- 
tions of H2 which are endothermic or possess activation 
barriers occur faster as the internal excitation of H2 in- 
creases. In some cases the internal energy of excited H2 
is effectively used to overcome the reaction barrier and 
the rate constant approaches the collision limit, whereas 
in other cases an activation barrier persists even if the 
internal energy of excited H2 makes the reaction to be 
exothermic. The rate constant enhancement is therefore 
quite specific of each reaction and difficult to predict 5 . In 
particular, it is not simply given by the subtraction of the 
internal energy of excited H2 from the endothcrmicity or 
activation barrier of the reaction with ground state H2 , as 
has b een assumed previously (e.g. iSternberg fc Dalgarnol 
1995). The rate constant enhancement when H2 goes 

5 We also note that for fast exothermic reactions (without bar- 
riers), the internal excitation of a reactant can actually inhibit 
the reactivity, as ob served experimentall y for ion-molecule (e.g . 
IGerlich & Roari989l ) and neutral-neutral l|01khov fc Smith! [20071) 
reactions. 



from the ground to an excited vibrational state gets 
larger at low temperatures and can reach many orders 
of magnitude. For instance, at 200 K it is somewhat less 
than a factor of ten for the reaction with CN, whereas 
for H2 + C + it is as high as ten orders of magnitude. 

The state-specific rate constants for reactions of H2 are 
obviously useless to study a gas in which the H 2 states 
are thermally populated, in which case the thermal rate 
constants are to be used, but become very useful for the 
study of regions where H2 levels have non thermal pop- 
ulations, which is found to occur quite often in the in- 
terstellar medium. In such regions some state-specific 
reactions of vibrationally excited H2 may become impor- 
tant to determine the chemical composition of the gas. 

Let us consider a region with a gas kinetic temperature 
T where H2 levels are not thermally populated, so that 
the fractional abundance of vibrationally excited H2 is /* 
= n(H2 v > 0)/n(H2 v = 0), and the generic reactions 

H 2 (v = 0) + X fe -^ } products (6) 

H 2 (v > 0) + X fc 3? products (7) 

so that Afc(T) = k*(T)/k°(T) is the rate constant 
enhancement at a temperature T. For reaction ([7]) to 
become important it has to be competitive with respect 
to reaction (j6]), i.e. f *Ak(T) > 1, and also with respect 
to other reactions which consume the reactant X. 
According to this, reactions such as H2 + CN, which has 
a moderate rate enhancement (Ak < 10), need a very 
high abundance of vibrationally excited H2 (/* > 0.1) to 
have some impact on the chemistry. On the other hand, 
reactions with a very high rate constant enhancement, 
such as that of H2 and C + for which Ak ~ 10 10 at 
200 K, just need a moderate abundance of vibrationally 
excited H2 (/* > 10~ 10 ). This is however a condition 
necessary but not sufficient to have some impact on the 
chemistry as the reaction rate of, in this case, H2(v > 0) 
+ C + needs also to be large enough compared to other 
reactions which consume C + . 

2.2. Implications for the interstellar medium 

We may state that the importance of the chemistry 
of vibrationally excited H2 is limited to reactions with 
a large rate constant enhancement and to astronomical 
regions with a high abundance of vibrationally excited 
H2. In order to get a quantitative estimate of the im- 
pact of these reactions on the chemistry of the interstellar 
medium we have run single point time dependent chem- 
ical models covering a wide range of the space of param- 
eters: fraction of vibrationally excited H2 (/*), density 
of hydrogen nuclei («h), far ultraviolet (FUV) radiation 
field strengt h (y, relative to the Draine interstellar ra- 
diation field: lDrainelll978() . and gas kinetic temperature 
(T). 

We consider 127 different species composed of the el- 
ements H, He, C, N, O, and S, whose ab undances are 
assumed to be solar (jAsplund et al.l [2009h with a 25 % 
of depletion for elements heavier than He, and assume 
that initially all the elements are in atomic form, ei- 
ther neutral or ionized. The chemical network com- 
prises 2307 gas phase reaction s whose rate constants 
have been mostly taken from [Agundez fc Cernicharol 
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Fig. 2. — Steady state abundances relative to total H nuclei for CH + , CtLJ" , CHg , CH, OH, and H2O as a function of njj/x f° r various 

= 10 3 cm" 3 (so that the FUV field strength \ has been 



fractional abundances /* of H2 (v 
varied), A v = 0.5, and T = 100 K 



1). All the models have been run with njj 



(2006) and from the UMIST database for astrochem- 
istry 6 (|Woodall et al.l [2007h . some of them revised ac- 
cording to the most recent literature on chemical kinet- 
ics. The p hotodissociation r ates of H2 and CO have been 
taken from lLee et al.l ([1996). The adopted visual extinc- 
tion is 0.5 mag and the cosmic-ray ionization rate per 
H2 molecule is 5xl0~ 17 s _1 . We do not consider grain 
surface reactions except for the formation of H2, which is 



sj, the canonical value in typical diffuse clouds (jJural 

This quite standard chemical network is completed by 
adding the reactions of vibrationally excited H2 detailed 
in Table [T] Thus, H2 (v = 1) is included as a new 
chemical species with a fractional abundance /*, rela- 
tive to ground state H 2 , which is kept fixed at four dif- 
ferent values (0, 10~ 6 , 10 -5 , and 10 ), independently 
of the mechanism responsible of the excitation of H2. 
The amount of H2 (v > 1) is important for reactions ^ 
and ([3]) so that we also consider H2 in the v — 2 and 
v = 3 states with populations which are 1/3 and 1/9 
respectively of that in the v = 1 state, values typically 
attained in interstellar clouds wit h an important fraction 
of vibrationally excited H2 (e.g. Ivan Dishoeck fc Black! 
ll986tlMever et aLl r2001). In these exploratory models we 
do not distinguish between the different reactivity of H2 
in the various rotational levels of the ground vibrational 
state. 

The abundance of various species are found to be sensi- 
tive to the fraction of vibrationally excited H2. In Fig. [5] 
we plot the steady state abundance of some of these 
species for various fractional abundances of Ha (t> = 1) 
as a function of the n^/x ratio. The steady state abun- 
dances are found to depend on this latter ratio rather 
than on either «h or x- For the results shown in Fig. [2] a 
density of hydrogen nuclei of 10 3 cm -3 has been adopted, 
so that the FUV field strength has been varied. Reac- 
tions of vibrationally excited H2 have a larger impact 

6 See http://www.udfa.net/ 



on the chemistry at low temperatures, since as discussed 
previously the rate constant enhancement Afc of the var- 
ious reactions is larger. The results shown in Fig. [^cor- 
respond to a relatively low temperature of 100 K. 

The molecule whose abundance is more affected by 
the reactions of vibrationally excited H2 is the methyli- 
dyne cation (CH + ), which for fractional abundances of 
H2(w = 1) above 10~ 7 — 10~ 6 is readily formed by the 
reaction between H2(v = 1) and C + . The steady state 
abundance of CH + approximately scales with the frac- 
tion of H2(w = 1) and reaches its maximum abundance 
for nn/x ratios in the range 1 — 100. By running mod- 
els at different kinetic temperatures we find that above 
400-500 K the reaction between H 2 (w = 0) and C+ 
starts to proceed and the abundance of CH + depends 
to a much lesser extent on the fraction of vibrationally 
excited H 2 . 

The abundances of CHj" and CH^" are strongly linked 
to that of CH + by the chain of exothermic reactions 



CH+ ^ CH+ % CH+ 



(8) 



which stops at CHjj", as this species does not react fast 
with H2. Therefore, in those conditions in which CH + 
reaches a relatively high abundance, CH^ and CH^" do 
so (see Fig. [2]). Some other species such as CH, OH, and 
H2O also experience a moderate abundance enhancement 
when CH + is formed. The radical CH is formed by the 
dissociative recombination with electrons of the cations 
CH^ and CH^, while OH and H2O are produced by a 
sequence of reactions that starts with the photodissocia- 
tion of CH + and production of H + followed by 



O 



0+ 5? OH+ 2? H 2 0+ 5? H 3 0+ ^ OH, H 2 (9) 



The reaction of vibrationally excited H2 and atomic 
oxygen has little effect on the abundance of OH for most 
of the conditions explored. It does contribute to OH 
formation just for fractional abundances of H2(w = 1) 
in excess of 10 -4 and for kinetic temperatures around 
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Fig. 3. — Abundances relative to the total number of H nuclei 
in C Oph for H2, H2 (v > 0), C+, and CH+ as a function of A v 
as given by the If^ model (dashed line corresponds to the CH+ 
abundance under the standard model). Also shown is the fraction 
of vibrationally excited H2 /* (referred to the left axis) and the 
gas kinetic temperature T (referred to the right axis). 

300 K. For much lower temperatures the activation bar- 
riers of the reactions H 2 (v > 0) + prevent to reach a 
high enough OH formation rate whereas for much larger 
temperatures the rate constant enhancements Afc are low 
and OH is mostly formed by the reaction of H 2 (v = 0) + 
O. As concerns the reactions of vibrationally excited H2 
with He + , OH, and CN, they have almost no impact on 
the chemistry as the rate constant enhancements are just 
moderate and require of a very large fraction of vibra- 
tionally excited H2 to compete with the corresponding 
reactions of H 2 (v = 0). 

In summary, from the five reactions of vibrationally 
excited H 2 considered here, only that with C + has a 
large impact on the abundances of some species, espe- 
cially CH + , in those astronomical regions which are not 
too warm (with kinetic temperatures below a few hun- 
dreds of degrees Kelvin) and which have a moderately 
high fraction of vibrationally excited H 2 (in excess of 
10~ 7 - 10~ 6 ). Reactions of H 2 (w > 0) with He+, O, OH, 
and CN are less likely to affect the chemistry of interstel- 
lar clouds unless large amounts of vibrationally excited 
H 2 are present. 

3. EFFECTS ON SELECTED ASTRONOMICAL SOURCES 

Vibrationally excited molecular hydrogen is commonly 
observed in PDRs, where the H 2 v > levels are excited 
by FUV fluorescence, and in shocked gas where the ex- 
citation is mainly collisional. Among these two types 
of regions the chemistry of vibrationally excited H 2 is 
expected to be less important in shocked regions since 
kinetic temperatures are usually very high and the H 2 
level populations a re not as far from thermalization as 
they are in PDRs ()Burtonlll992D . Hereafter we thus fo- 
cus on the case of PDRs to investigate the effects of the 
reactions of vibrationally excited H 2 onto the gas chem- 
ical composition. 

For that purpose we have utilized an updated version 
of the Meudon PDR code 7 , a phot ochemical model of a 
plane parallel and stationary PDR dLe Petit et alj|2005 
Goicoechea fc Le Bourlotj 120071 : iGonzalez Garcia et al.l 
20081) . to compute the chemical and physical structure 

7 See http://pdr.obspm.fr/PDRcode.html 



as a function of the cloud depth for various PDR-like 
regions. Penetration of FUV radiation strongly de- 
pends on dust grain properties through dust absorp- 
tion and scattering of FUV photons which, in addi- 
tion, determine the efficiency of the dominant heat- 
ing mechanisms of the gas (photoelectric heating). In 
this work we use standard dust pr operties, i.e. grain 
sizes follow a power-law distrib ution (M athis et al.lll977t 
iGoicoechea fc Le Bourlotl[2007D and explicitly determine 
the gas temperat ure gradient by solvin g the cloud ther- 
mal balance (see iLe Petit et al.l 120061 ) . We have used 
the chemical network described in the previous section 
and the reactions of excited H 2 shown in Table [1] For 
the reaction between H 2 and C + we adopt the rotational 
state-specific rate constant expressions (No 2 in Table [lj 
for H 2 (j <7) and assume the Langevin value (No 3 in 
Table [IJ for H 2 (j > 7 or v > 0), states for which the 
reaction becomes exothermic. Under thermal conditions 
this approach results in a value equivalent to the thermal 
rate constant (No 1 in Table [lj, at worst a factor of ^2 
larger in the high temperature range (>I000 K), where 
the rate constant starts to be dominated by H 2 v > 
states. 

We consider (de-)e xcitation of H 2 levels b y collisions 
with H 2 , H, and He (|Le Bourlot et al.lll999D . Collisions 
with H + are a lso included based on the recent work of 
Huestis (2008), although they are found to be not im- 
portant. Concerning collisions with electrons, the rates 
of vibrational de-excitation of H 2 v = I — > are typi- 
cally a few 10 ~ 10 cm 3 s" 1 (fro m the cross sections rec- 
ommended bv lYoon et a l. 2008), which implies a critical 
density of electrons of ~ I0 3 cm -3 . In low density re- 
gions, such as diffuse clouds, the electron density is much 
lower and electron-impact vibrational excitation is neg- 
ligible, although in denser regions with a high ionization 
degree, such as planetary nebulae, electrons could play 
a role in producing non-thermal H 2 vibrational popula- 
tions. The implementation of (de-)excitation of H 2 by 
collisions with electrons is however not straightforward 
as only a few sta te-to-state cross sections are available 
HYoon et al.ll2008l) . 

We are mainly interested in investigating how impor- 
tant is the effect of the chemistry of vibrationally excited 
H 2 in the various regions. We thus have systematically 
run for the different sources a model in which we adopt 
the thermal rate constant expressions for all reactions 
of H 2 (standard model) and another one including the 
state-specific rate constants for the reactions of excited 
H 2 detailed in Table Q] (if 2 model). 

3.1. Diffuse clouds: ( Oph and ED 34078 

We first focus on a well known diffuse cloud such as 
£ Op h, whose parameters (s e e Tab le [2t have been taken 
from Ivan Dishoeck fc Blackl (fl986l) . In Fig. [3] we show 
as a function of Ay the abundance of some species while 
the calculated column densities are given in Table [5J The 
largest effects of the chemistry of vibrationally excited 
H 2 are seen at the edges of the cloud (Ay < 0.1) where 
there is an important fraction of H 2 (v > 0), which gives 
rise to an increase in the abundance of CH + . Note that 
the sharp increase of /* (the ratio of vibrationally ex- 
cited to ground state H 2 ) up to ~ 10~ 3 at the very outer 
edges (Ay < 0.001) is produced by the decrease in the 
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TABLE 2 

Calculated column densities of selected species in the studied astronomical sources 



rin ( cm 3 ) 

^total 



HD 34078 "hot PDR" 
F3 



"17 



rion Bar 



NOC 7027 



W 
3xl0 4 || 
9.0 
20-700 



8x10° 
5xl0 3b || 
2.5 
200-800 



HD 176386" 
3xl0 3 
125 c || 

0.9 
95-120 



600 
4 || 4 

0.87 
40-50 



10* 
10 4 || 

1.0 
200-300 





standard - 




standard - 




standard 


— > 








standard 


-» 


H* 




standard — ¥ H% 


H 


1.2xl0 au - 


-> ~ 


1.6xl0 2i - 




1.4xl0 21 


— »• 








2.1X10 21 








l.lxlO 2 " - 


j. ~ 


H 2 


7.6 xlO 20 - 




1.6X10 20 - 




7.7X10 21 


— > 








1.3X10 21 


— > 






8.8X10 20 - 


-¥ 7.9xl0 20 


H 2 (> > 0) 


3.7xl0 13 - 


* ~ 


4.2X10 15 - 




2.2X10 16 










4.9xl0 17 


— > 


4.2x 


10 17 


1.8X10 14 - 


T> ~ 


C+ 


5.6xl0 17 - 


•> ~ 


2.5xl0 17 - 




7.5xl0 17 


-> 








4.6xl0 17 


— > 


3.4x 


10 17 


3.6xl0 17 - 


-> ~ 


CH+ 


5.2xl0 10 - 


-5- 9.3xl0 10 


1.4xl0 10 - 


-> 1.6xl0 12 


4.1X10 11 


— ► 


l.l 


X 


10 13 


5.7xl0 13 


— > 


9.6x 


10 13 


5.1xl0 10 - 


-»■ 1.9X10 11 


CH+ 

chJ 


1.3X10 11 - 


•> ~ 


3.3xl0 10 - 


-> 6.5X10 11 


4.6X10 11 


— > 


6.6 


X 


10 12 


1.8xl0 13 


— > 


3.5x 


10 13 


1.3X10 11 - 


-¥ 1.9X10 11 


2.3X10 11 - 




3.8xl0 10 - 


-> 6.6X10 11 


1.9xl0 12 


— > 


2.1 


X 


10 13 


5.0xl0 13 


— > 


2.0x 


10 14 


4.7xl0 n - 


-> 2.0X10 11 


CH 


1.7xl0 13 - 


-> ~ 


2.6xl0 10 - 


-> 4.4X10 11 


4.3xl0 13 


— > 


8.6 


X 


10 13 


7.7xl0 13 


— > 


8.4x 


10 13 


3.2xl0 12 - 


-> 4.8xl0 12 


OH 


l.lxlO 12 - 


j. ~ 


3.7xl0 10 - 


4 3.1X10 11 


4.6xl0 13 


— > 








2.5xl0 15 


— > 






l.lxlO 11 - 


4 4.4X10 11 


H2O 


l.OxlO 11 - 


•> ~ 


2.5xl0 9 - 


-> 2.3xl0 10 


4.4xl0 15 


— > 








l.OxlO 15 


— > 






l.lxlO 10 - 


4 4.7xl0 10 


CO 


8.0xl0 12 - 


* ~ 


8.2xl0 9 - 


-> 1.6X10 11 


1.4xl0 18 


— y 








2.2xl0 18 


— > 






2.0xl0 12 - 


3.0xl0 12 


CN 


4.4xl0 n - 


■> ~ 


1.0x10 s - 


-> 2.4xl0 9 


2.1xl0 13 


— > 


2.6 


X 


10 13 


3.2xl0 15 


— > 






1.7xl0 10 - 


4 4.4xl0 10 


HCN 


5.4xl0 8 - 


-> ~ 


1.1X10 6 - 


-> 2.4xl0 7 


1.5xl0 13 


-> 








7.6 xlO 15 


— > 






9.2xl0 7 - 


-> 1.7x10 s 



Note. — The adopted elemental abundances for C, N, O, and S are respectively 3.51x10 , 7.5x10 , 6.23x10 , and 9.9xl0~ 6 in C Oph; 
1.32X10 -4 , 7.5xl0~ 5 , 3.19X10 -4 , and 1.85X10 -5 in the "hot PDR" of HD 34078; 1.4xl0~ 4 , l.OxlO -4 , 3.0xl0 -4 , and 2.8X10" 5 in Orion Bar; 
1.26xl0~ 3 , 1.91X10 -4 , 5.5xl0~ 4 , and 7.94xl0" 6 in NGC 7027; and 2.02xl0~ 4 , 5.07xl0" 5 , 3.67xl0" 4 , and 9.89xl0~ 6 in HD 176386. 
a The notation x \ | y refers to the FUV field x and y at each side of the cloud. 

b FUV field at 0.0169 pc of a star with T» = 198000 K and R, = 0.075 R radiating as a blackbody. 
c FUV field at 0.04 pc of a star with X* = 11300 K and i?, = 2.55 Rq radiating as a blackbody. 
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Fig. 4.— Same as in Fig.[3]but for the "hot PDR" of HD 34078. 



abundance of H2 (v = 0), relative to total hydrogen, 
while that of H 2 (v > 0) keeps roughly constant. In this 
regime large /*A fc(T) ratios are attained (see discussion 
at the end of Sec. 12.11) and reactions of vibrationally ex- 
cited H2 increase their importance with respect to those 
of ground state H2, but not necessarily with respect to 
other processes (e.g. photodissociations). That is why 
the jump of /* at the outer edge of £ Oph is not followed 
by the abundance of CH + (see Fig. [3J. 

In any case, as the fraction of H2 (v > 0) remains rela- 
tively low (/* < 10 -7 ) throughout most of the cloud, 
the total column density of CH + is not dramatically 
affected. It increases by about a factor of 2 (see Ta- 
ble [2]) but remains well below the observed value of 
2.9 xlO 13 cm" 2 (see Ivan Dishoeck fc Black1ll986D . This 
relatively low fraction of vibrationally excited H2 is 
in agreement with ast ronomical observations of £ Oph 
(|Federman et a l. 1995), and implies that for quiescent 



diffuse clouds the chemistry of vibrationally excited H2 
does not have dramatic effects on the global chemical 
composition of the cloud. 

Large fractions of vibrationally excited H2 are not usu- 
ally observed in diff use clouds, except for a few lines 
of sight: HP 37903 (IMever et al.ll200l and HD 34078 
( Boi sse et all l200l |2009j), although its presence seems 
to be related to a dense component exposed to a high 
FUV field arising from nearby hot stars. In Fig. fj] we 
show the results obtained for the dense component o f 
HD 34078, called the "hot PDR" bv lBoisse etall (|2005D . 
whose parameters have been taken from these authors 
(see Table [2]). In this case, due to the large FUV field, a 
high fraction of H2 (v > 0) is maintained throughout the 
cloud. Thus, the reaction between H2 (v > 0) and C + 
readily occurs affecting the abundance of various species 
(see column densities in Table H]), most notably CH + 
whose abundance is increased by two orders of magnitude 
with respect to that obtained in the standard model, yet 
remaining below t he observed value (6.5xl0 13 cm" 2 ; see 
IBoisse et al.ll2005l) . We note that the differences between 
the column densities calculated in our standard model 
and those obtained bv lBoisse et al.l ([20 05^ are most likely 
due to differences in the chemical network utilized. For 
instance, the larger column density of CH + obtained by 
them is to a large extent due to the larger rate constant 
they adopt for the radiative association between C + and 
H2 : 1.7x 10~ 15 c m 3 s" 1 (upper limit derived for para H2 
bv lCxerlicU[l99l versus the more conservative value of 
4xlQ- 16 (r/30 0)~°- 2 cm 3 s" 1 ado pted by us (based on a 
calculation bv lHerbst et al.lll985l ). 

We thus find that the reaction of H2 (v > 0) and C + is 
not able by itself to explain the large CH + column den- 
sities observed in diffuse and translucent clouds (> 10 13 
cm" 2 ). Nowadays, the most popular theories on the CH + 
formation in diffuse clouds invoke the existence of tcm- 
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Fig. 5. — Same as in Fig.[3]but for the Orion Bar. Fig. 6. — Same as in Fig.[3]but for NGC 7027. 



porar y events, such as shocks (iPineau des Forets et al.l 
119861) or turbulence dissipation (jGodard et aiT l2009). 
which increase locally the gas kinetic temperature and 
allow rea ction (HI) to proceed. An other possibility is that 
given by iFederman et all (|T996), who suggest that the 
ion-neutral drift produced by Alfven waves could result 
in a non thermal rate enhancement of ion-neutral reac- 
tions, allowing CH + to be effectively formed through re- 
action (jTJ) . An enhanced FUV field has also bee n pro- 
posed as the ultimate cause of CH + formation (jSnowl 
1993), hypothesis that is supported by the fact that 
in diffuse clouds the column density of CH + is well 
correlated with rotation ally excited H2 (J » 3 — 5) 
([Lambert fc Danks 198(|), levels which are to a large ex- 
tent pumped by FUV fluorescence in diffuse clouds. An 
intense FUV field will indeed bring H2 molecules to vi- 
brational excited states, which in turn will react fast with 
C + to form CH + , but when the various magnitudes are 
evaluated through a PDR model under a plausible phys- 
ical scenario for a diffuse cloud it is found that the abun- 
dance of H2 (v > 0) is not high enough to form CH + with 
an abundance close to the the typical observed values, as 
shown in the case of C Oph. 

3.2. Dense interstellar PDRs: Orion Bar 

Dense and highly FUV-illuminated clouds are regions 
where the chemistry of vibrationally excited H2 is ex- 
pected to play an important role as the fraction of H2 
(v > 0) is usually large. Here we consider the case of one 
of the most archetypical dense PDRs which is the Orion 
Bar, whose parameters ( see T able [2J have been taken 
from Ivan der Wiel et all (pOOl . In this case the combi- 
nation of a high FUV radiation field and a high density 
( n ~B./x ~ 3) results in a fraction of H2 (v > 0) in excess of 
10 -6 in the region Ay — 0.5 — 1.5, where hydrogen is still 
mostly molecular and the gas kinetic temperature is just 
of a few hundreds of degrees Kelvin (see Fig. [5]). Thus, 
the reaction H 2 (v > 0) + C + has indeed a great impact 
on the abundance of CH + which increases by almost 2 
orders of magnitude with respect to that obtained in the 
standard model. 

The existence of an important amount of vibrationally 
excited H2 in th e Orion Bar is supported by astronomi- 
cal observations (Ivan der Werf et~atlll996D . On the other 
hand, CH + has not yet been observed but its pure ro- 
tational lines lie in the wavelength range of the Herschel 



Space Observatory. If, as indicated by our model, the 
formation of CH+ occurs by the reaction of H 2 (v > 0) 
and C + , then the spatial distribution of CH + emission 
should be well correlat ed with that of the v = 1— >0 S(l) 
line of H2 observed by Ivan der Werf et all (fl996h 

3.3. Planetary nebulae: NGC 7027 

Planetary nebulae are also sources where the chemistry 
of vibrationally excited H2 can be important since the cir- 
cumstellar gas is normally exposed to a very high FUV 
flux which emanates from the hot central white dwarf. A 
prototype of this kind of objects is NGC 7027, which har- 
bors one of the hottest known stars (T* ~ 200,000 K). To 
model this source we have adopte d the physical scenario 
utilized bv lHasegawa et all ()2000D focusing on the dense 
region (see Table [5]). As shown in Fig. [5] a large fraction 
of H2 (v > 0) is present in the circumstellar gas of this 
source (/* ~ 10~ 4 — 10~ 3 ), and CH+ is readily formed 
by the reaction of H2 (v > 0) and C + in the regions lo- 
cated close to the star, where C + reaches its maximum 
abundance. In these regions the reaction H2 (v = 0) + 
C + is also an important source of CH + due to the high 
temperatures attained. 

From the observational point of view the existence 
of vibrationally excited H2 in N GC 7027 is well es- 
tablis hed from ISO observations (|Bernard Salas et"al"1 
l200l . as is the presence of CH + with a column density 
of 0.8-2.5xl0 14 (ICernicharo et al.lll997l) . in agreement 
with the value we get. 

3.4. Circumstellar disks: HD 176386 

Finally, we treat also the case of circumstellar disks 
surrounding young stars of T Tauri and Herbig Ac/Be 
type, where the outer disk surface located at just a few 
AU of the central star is exposed to the strong stel- 
lar UV flux. Hot molecular hydrogen both rotationally 
and vibrationally e xcited, has been observed in various 
sources of this type (jMartin-Zaidi et al.ll2008l ). To inves- 
tigate the importance of the chemistry of vibrationally 
excited H2 in such environments we take the parame- 
ters derived for the warm H2 emitting gas observed by 
iMartin-Zai'di et al.l (|2008l ) toward the B9 star HD 176836 
(see Table [2]). As shown in Fig. [7] the fraction of vibra- 
tionally excited H2 stays around 10~ 7 — 10~ 6 and so the 
abundance of CH + is increased by almost one order of 
magnitude when reactions with excited H2 are included, 
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Fig. 7.— Same as in Fig.[3]but for HD 176386. 

although in the case of HD 176386 the total column den- 
sity (see Table 0) remains probably too low for being 
detected. Some other species experience moderate abun- 
dance enhancements as a consequence of the activation 
of the reaction between H2 (v > 0) and C + . 

4. SUMMARY 

In this paper we have investigated the impact of vari- 
ous chemical reactions involving vibrationally excited H2 
which have been studied either experimentally or theo- 
retically (concretely those with C + , He + , O, OH, and 
CN) on the chemical composition of various astronomical 
regions. Among the reactions considered here that of H2 
(v > 0) with C + stands out as the most important one. 
This reaction becomes the main formation route of the 
reactive cation CH + and controls the abundance of some 
other related species in moderately warm astronomical 
environments with fractions of vibrationally excited H2 
in excess of some 10 -6 , conditions which are found to 
occur in some types of PDR-like regions. 

The importance of the chemistry of vibrationally ex- 
cited H2 has been investigated through PDR models in 
the diffuse clouds C Oph and HD 34078, the Orion Bar, 
the carbon-rich protoplanetary nebula NGC 7027, and 
the circumstellar disk HD 176386. The study of £ Oph 
indicates that the fraction of H2 (v > 0) present (< 
10~ 7 ) is insufficient to form CH + with the abundance 
typically observed in diffuse and translucent clouds. On 



the other hand, in dense and highly FUV illuminated 
clouds, such as the hot PDR of HD 34078, the Orion 
Bar and NGC 7027, reactions of vibrationally excited 
become crucial to determine the global chemical com- 
position. In particular, the reaction H2 (v > 0) + C + 
becomes the dominant synthetic pathway to CH + and 
make the related species CH^ and CH^ to reach similar 
or slightly lower abundances. These two cations however 
have not yet been observed in any astronomical region — 
the spectroscopic data of CH^ is poorly known (see e.g. 
IPolehampton et all 1200 71 ) while CH;j~ lacks a permanent 
dipole moment. 

Other reactions of vibrationally excited H2 not yet 
studied either experimentally or theoretically, and thus 
not included here, could become important in regulating 
the chemical balance of interstellar clouds. For instance, 
the reactions H 2 + C -)• CH + H, H 2 + S -> SH + H, 
and H 2 + S+ -)• SH+ + H are endothermic by 0.99 eV 
(~ 11500 K), 0.83 eV (~ 9600 K), and 0.86 eV (~ 9860 
K) respectively, and their rate constants have an expo- 
nential term of about the value of the endothermicity. 
This may indicate that, analogously to reaction ([1]), the 
internal energy of excited H2 would be effectively used to 
diminish or overcome the activation barrier so that for 
H2 (v > 1), all these reactions being exothermic, they 
would proceed at about the collision limit. Experimental 
or theoretical information on these or other potentially 
important reactions of vibrationally excited H2 would be 
of great interest. 

In any case we encourage astrochemists to include the 
reactions of excited H2 for which kinetic data is available 
(see Table [l| when modelling the chemistry of regions 
with an important fraction of vibrationally excited H2. 
At least the reaction between H2 (v > 0) and C + clearly 
has an important impact on the global chemical compo- 
sition of such regions. 



We thank the referee for constructive comments on 
this article. M.A. is supported by a Marie Curie Intra- 
European Individual Fellowship within the European 
Community 7th Framework Programme under grant 
agreement n° 235753. J.R.G. was supported by a Ramon 
y Cajal research contract from the Spanish MICINN and 
co-financed by the European Social Fund. 



REFERENCES 



Agundez, M. & Cernicharo, J. 2006, ApJ, 650, 374 
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, 

ARA&A, 47, 481 
Bernard Salas, J., Pottasch, S. R., Beintema, D. A., & Wesselius, 

P. R. 2001, A&A, 367, 949 
Boisse, P., Le Petit, F., Rollinde, E., et al. 2005, A&A, 429, 509 
Boisse, P., Rollinde, E., Hily-Blant, P., et al. 2009, A&A, 501, 221 
Burton, M. G. 1992, Aust. J. Phys., 45, 463 
Cernicharo, J., Liu, X.-W., Gonzalez-Alfonso, E., et al. 1997, 

ApJ, 483, L65 

Chase, M. W. 1998, NIST-JANAF Thermochemical Tables, 4th 

ed., J. Phys. Chem. Ref. Data, Monograph n 9 
Dabrowski, I. 1984, Can. J. Phys., 62, 1639 
Draine, B. T. 1978, ApJS, 36, 595 

Federman, S. R., Cardelli, J. A., van Dishoeck, E. F., Lambert, 

D. L., & Black, J. H. 1995, ApJ, 445, 325 
Federman, S. R., Rawlings, J. M. C, Taylor, S. D., & Williams, 

D. A. 1996, MNRAS, 279, L41 
Freeman, A. & Williams, D. A. 1982, Ap&SS, 83, 417 



Garrod, R. T., Rawlings, J. M. C, & Williams, D. A. 2003, 

Ap&SS, 286, 487 
Gcrlich, D., Disch, R., & Scherbarth, S. 1987, J. Chem. Phys., 87, 

350 

Gerlich, D. & Rox, T. 1989, Zeit. Phys. D: Atoms, Molecules and 

Clusters, 13, 259 
Gerlich, D. 1994, in AIP Conf. Proc. 312, Molecules and Grains 

in Space, ed. I. Nenner (New York: AIP), 489 
Godard, B., Falgarone, E., & Pineau des Forets, G. 2009, A&A, 

495, 847 

Goicoechea, J. R. & Le Bourlot, J. 2007, A&A, 467, 1 
Gonzalez Garcia, M., Le Bourlot, J., Le Petit, F., & Roueff, E. 

2008, A&A, 485, 127 
Hasegawa, T., Volk, K., & Kwok, S. 2000, ApJ, 532, 994 
Herbst, E. 1985, ApJ, 291, 226 

Hierl, P. M., Morris, R. A., & Viggiano, A. A. 1997, 

J. Chem. Phys., 106, 10145 
Hollenbach, D. & McKee, C. F. 1989, ApJ, 342, 306 
Huestis, D. L. 2008, Planet. Space Sci., 56, 1733 



The chemistry of vibrationally excited H 2 



9 



Jones, M. E., Barlow, S. E., Ellison, G. B., & Ferguson, E. E. 

1986, Chem. Phys. Lett., 130, 218 
Jura, M. 1975, ApJ, 197, 575 

Lambert, D. L. & Danks, A. C. 1986, ApJ, 303, 401 
Le Bourlot, J. 1991, A&A, 242, 235 

Le Bourlot, J., Pineau des Forets, G., & Flower, D. R. 1999, 

MNRAS, 305, 802 
Le Petit, F., Nehme, C, Le Bourlot, J., & Roueff, E. 2006, ApJS, 

164, 506 

Lee, H.-H., Herbst, E., Pineau des Forets, G., Roueff, E., & Le 

Bourlot, J. 1996, A&A, 311, 690 
Martin-Za'idi, C, Deleuil, M., Le Bourlot, J., ct al. 2008, A&A, 

484, 225 

Mathis, J. S., Rumpl, W., Nordsieck, K. H. 1977, ApJ, 217, 425 
Meyer, D. M., Lauroesch, J. T., Sofia, U. J., Draine, B. T., & 

Bertoldi, F. 2001, ApJ, 553, L59 
Olkhov, R. V. & Smith, I. W. M. 2007, J. Chem. Phys., 126, 

134314 

Park, I. H., Wakelam, V., & Herbst, E. 2006, A&A, 449, 631 
Pineau des Forets, G., Flower, D. R., Hartquist, T. W., & 

Dalgarno, A. 1986, MNRAS, 220, 801 
Polchampton, E. T., Baluteau, J. -P., Swinyard, B. M., ct al. 2007, 

MNRAS, 377, 1122 



Schauer, M. M., Jefferts, S. R., Barlow, S. E., & Dunn, G. H. 

1989, J. Chem. Phys., 91, 4593 
Snow, T. P. 1993, ApJ, 402, L73 

Stecher, T. P. & Williams, D. A. 1972, ApJ, 177, L141 
Sternberg, A. & Dalgarno, A. 1995, ApJS, 99, 565 
Sultanov, R. A. & Balakrishnan, N. 2005, ApJ, 629, 305 
Tielens, A. G. G. M. & Hollenbach, D. 1985, ApJ, 291, 722 
Truong, T. N. 1995, J. Chem. Phys., 102, 5335 
van der Werf, P. P., Stutzki, J., Sternberg, A., & Krabbe, A. 

1996, A&A, 313, 633 
van der Wiel, M. H. D., van dcr Tak, F. F. S., Ossenkopf, V., et 

al. 2009, A&A, 498, 161 
van Dishoeck, E. F. & Black, J. H. 1986, ApJS, 62, 109 
Wagner, A. F. & Graff, M. M. 1987, ApJ, 317, 423 
Woodall, J., Agundez, M., Markwick-Kcmpcr, A. J., & Millar, T. 

J. 2007, A&A, 466, 1197 
Yoon, J.-S., Song, M.-Y., Han, J.-M., ct al. 2008, 

J. Phys. Chem. Ref. Data, 37, 913 
Zellner, R. & Steinert, W. 1981, Chem. Phys. Lett., 81, 568 
Zhu, W., Zhang, J. Z. H., Zhang, Y. C, ct al. 1998, 

J. Chem. Phys., 108, 3509 



